Equations of Stellar Structure 


Stellar structure and evolution can be calculated via a series of 
differential equations involving mass, pressure, temperature, and 
density. For simplicity, we will assume spherical symmetry. 


First, 





let’s define the global variables 


r = the radial coordinate of the star 


7 = Total Mass of Star 

r = Total Luminosity of Star 

R= the Radius of the star 

t = the independent time variable 





Obviously 


Note: 


X+Y+Z=1 


I’ve written the above parameters in terms the radial co- 


ordinate of the star. However, for many reasons, it is better to 
use the Lagrangian formulation, and use M, the enclosed mass, 
and the independent variable. ‘The dependent variables are then 


r(M), £(M), P(M), etc. 





In addition, we will also define some local variables: 


Ky 


= the luminosity produced per unit mass of material 
= the opacity (photon cross section) per unit mass of material. 


Conservation of Mass 


For a spherically symmetric star, the mass interior to some radius 
Ris 


R 
M(r) = f Arr? p(r) dr (2.1.1) 
0 
Written in terms of the differential, this is 
d 
mM = Arr? p(r) (22) 


However, over its lifetime, a star’s radius will change by many 
orders of magnitude, while its mass will remain relatively con- 
stant. Moreover, the amount of nuclear reactions occuring inside 
a star depends on the density and temperature, not where it is 
in the star. ‘Thus, as stated above, a better and more natural 
way to treat stellar structure is to write radius as a function of 
mass, 1.e., 

dr 1 

dM  4rr?p 


This is the Lagrangian form of the equation (rather than the 
Eulerian form). All the equations of stellar structure will be 
expressed in the Lagrangian form, and most of the parameters 
will be expressed in per unit mass, rather than per unit size or 
volume. 


(2.1.3) 


Conservation of Momentum 


If a star is in Hydrostatic Equilibrium, then the force of gravity 
at any point must be balanced by the pressure. Consider a shell 
at radius r: the mass per unit area in the shell is pdr, and the 
weight (per unit area) is —godr. This weight must cancel the 
pressure difference felt from one side of the shell to the other, 
Pinner ~ Pouter = IÈ dr. Thus 


dP 


can gree 2.2.1 
a gp (2.2.1) 


Since the force of gravity, g is just GM/r?, then 





dP GM(r) 
= =P (2.2.2) 


So if we multiply by the conservation of mass equation, we get 
the Lagrangian form of hydrostatic equilibrium 











dP dr GM 1 
dp dM ( r? p) F. ae 
or 
dP GM 
I Ar 12 


Note, however, that stars need not be in hydrostatic equilibrium. 
Some stars pulsate! In these cases, the pressure at any moment 
may not equal the force of gravity, causing the shell to undergo 
an acceleration. 


Again, consider a shell at radius r. The force per unit area acting 
upon this shell is simply the difference in pressure from one side 
of the shell to the other 





dP 


while the force per unit area due to gravity is 


GM(r 
jo=gopar = oe), dr (2.2.6) 





These two forces define the acceleration of the shell 


fe+f,=- (Zar) + (Hoar) = pdr (=) (2.2.7) 





r2 


Again, this should be translated to the Lagrangian form by sub- 
stituting dM for dr using the mass conservation equation 


dP GM dM d?r 
(2.2.8 


(aii) =~ (ier) (ae) (Gee) 8 


Note that if the acceleration is zero, the equation reduces to the 
simple hydrostatic equation. 








Conservation of Energy 


Consider the net energy per second passing outward through 
a shell at radius r. If no energy is created in the shell, then 
the amount of energy in equals the amount of energy out, and 
d£L/dr = 0. However, if additional energy is created or absorbed 
within the shell, then d£/dr will be non-zero. Let’s define € as 
the energy released per second by a unit mass of matter. Then 





dL 
dE = 4rr?pe (2.3.1) 
or, in the Lagrangian form 
dL 


Note that e has three components. 

1) €,, the total energy created by nuclear reactions 

2) €,, the energy input into neutrinos, and 

3) €g, the energy produced or lost by gravitational expansion 
or contraction. 


Thus, 


dL 

dM 
In general, the contribution from nuclear reactions will always be 
positive, while the energy in neutrinos will always be lost from the 
system. (That’s why there’s a negative sign before e. But note 
that supernovae are the exception to this rule.) The gravitational 
energy can be either positive or negative, and simply reflects the 
heat gained or lost in the shell by PdV work. Thus, 


dq 
Eg = TE 


= €n — €y + €g (2.3.3) 





which, from thermodynamics (1.18), is 


_dq__ aT bP 
dt Td pdt 


Eg = 


(2.3.4) 


This can sometimes be written in a more useful form by recalling 
that 





Pô 
TONE 1.1 
Vad = a (1.19) 
eae ldr  NVsadP 
ee Ty meee 23: 
a (4 dt P =) eee) 


Thermal Structure 


The thermal structure of the star will depend on its method for 
transporting energy from its interior to its surface. For conve- 
nience, we will express the temperature stratification of the star 
in terms of the temperature gradient 


dln T 
= 241 
y (T) ( ) 


This temperature structure is then simply 





dT dP dT dP T 
na ae Ce) 


If the star is in approximate hydrostatic equilibrium, this simpli- 


fies to 
dT GMpT 


dr r2P 


or, in Lagrangian terms 


y (2.4.3) 





dT GMT 


dM‘ Anr4P ae) 


